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ABSTRACT 

Aims. To render cosmic shear an astronomical tool of high precision, it is essential to eliminate systematic effects upon its signal, one 
of the most significant ones being the intrinsic alignment of galaxies. The alignment in tidal fields that are created by the surrounding 
matter structure induces correlations between the intrinsic ellipticities of source galaxies, as well as correlations between the gravita- 
tional shear and the intrinsic ellipticity. While the former effect is restricted to physically close galaxy pairs and thus relatively easy 
to control, shear-ellipticity correlations occur for pairs at large separations. Because of the crudeness of current models of intrinsic 
alignment, we have developed a model-independent, purely geometrical method for removing the contamination of the cosmic shear 
signal by shear-ellipticity correlations. 

Methods. We remove the contributions to a tomographic cosmic shear signal that may be subject to contamination by shear-ellipticity 
correlations, making use of the characteristic dependence of these correlations on redshift. By introducing an appropriately chosen 
weight function to the lensing efficiency that nulls signals stemming from certain distances, new second-order measures of cosmic 
shear can be constructed that are free from intrinsic alignment. We present three approaches to determining such weight functions, 
optimized with respect to the amount of information the weighting preserves. After generalizing the construction of weight functions, 
the loss of information induced by this nulling technique and the subsequent degradation of constraints on cosmological parameters 
is quantified in a likelihood analysis. 

Results. For constructing optimal weight functions, good agreement is achieved between all approaches considered. In particular, a 
simplified analytical ansatz is shown to approximate the numerical results closely, significantly lowering computational efforts. For a 
survey divided into 20 redshift bins, we find that the area of credible regions increases by 20 % up to about 50 % after the application of 
nulling, depending on the cosmological parameters considered. We demonstrate that, due to the optimization of the weight functions, 
nearly all information is contained in a small subset of the new second-order measures. The use of a significantly smaller number of 
redshift bins than 20 for the nulling considerably degrades parameter constraints under conservative assumptions, emphasizing the 
need for detailed redshift information. 

Key words, cosmology: theory - gravitational lensing - large-scale structure of the Universe - cosmological parameters - methods: 
data analysis 



^ 1 . Introduction One of the major issues in rendering cosmic shear a tool of 

precision cosmology is the identification and removal of poten- 

Cosmic shear, the gravitational lensing of distant galaxies by tial systematic effects that could have an influence on the shear 

the intervening large-scale matter structure, provides a unique signal if measured to high accuracy. Apart from the intricacies of 

tool for cosmology, enabling cosmological parameters to be con- systematic errors in observations, one has to take approximations 

strained based on only a few physical assumptio ns. Detected made in the theory into account. Among these are the Born ap- 

for the first time at the turi i of the century by iBacon et al] proximation, the dropping of lens-lens coupHng and the replace- 

(2000), Kaiser et al. (2000j), Ivan Waerbeke et all (l2000|), and ment of the reduced shear by the shear itself, all fliree effect s 

[Wittman et al..(i200a) . cosmic shear observations have nowadays being a few per cent in magnitude jSchneider et alJll998ll2002l) . 

reached statistical errors on cosmological param eters compatible Physical effects to be considered are for instance source-lens and 

to Qiose of oth er estabHshed mefliods (see e.g. Jarvis et al. '2006'; source-redshift clustering, where in bo th cases the contaminat ion 



Hoekstra et al. 2006; Semboloni et al. 2006; Hetterscheidt et al. does not exceed the per cent level (see lSchneider et al.ll2002l) . 



2007; Beniamin et al. 2007; Fu et al. 2008). Lensing delivers in- a . n • j-c ci. ■ u 

J. . - , • • r„ T A potentially more serious modification of the cosmic shear 

formation complementary to CMB amsotropies, Type la super- ■ ,• . , , r - ■ i n- 

, , , ./ , — "^^^V — ^1 r.nnr^ \, signal IS caused by correlations of intrinsic galaxy ellipticities, 

novae, or galaxy redshift surveys (e.g. ISpergel et a l. 2007; Hu • • , ■ £c ^ r u w ^ j i . 

rrrrrr- 1 , , , . • , , • , , — \ — 2 mimicking the enect of shear correlations. Consider a correlator 

2002aD . thereby playing a crucial role in the development of pre- , ,,. u- u f »u u f n j 

' . . ^ , T • 1 1 ^ . of two galaxy ellipticities, which forms the base of all second- 

cision cosmology. In particular, cosmic shear surveys of the near , • . rr, / i \ j n- 

^ , ■, , , , • ■ order cosmic shear measures. The (complex) measured elliptic- 

future are currently considered to be among the most promising . ^, r • . • • n- s r i j 

, J, , , , . , Itt lAr^o^i lItt 7 J ity e IS the sum of the intrinsic ellipticity e of the galaxy and 

probes of dark energy and its parameters (see H ul2002bt Huteren i u ^ ^ c ^ ^ ■ 

r|;„„- . , ^ \r^r^r i iV-i^^^T^ thc gfavitational shcar Y. Writing thc coiTcIator for two galaxics 

12002; Albrecht et al. 2006; Peacock et al. 2006,) . ■ 5 ■ • u 

I and } yields 
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For the following discussion assume that galaxy / is situated 
closer to the observer, i.e. z, < zj- 

The first term on the righthand side is the desired cosmic 
shear signal. The second term might be non-vanishing if the two 
galaxies were subject to the same tidal forces during their for- 
mation or evolution, for instance if they had formed in a single 
dark-matter halo. This can only happen if the galaxies are phys- 
ically close, both on the sky and along the line of sight, so if 
Zj - Zi « 1. In the following we use the name 'intrinsic elliptic- 
ity correlation' for this effect. 

Various publi c ations (e.g. ICroft & Metzlei 120001: 
ICatelan et al.l 1200 U ICrittenden et al.l l2002t iJind |2002|) with 
both analytical and numerical approaches, e.g. via N-body 
simulations, have dealt with this possible contaminant of cosmic 
shear They agree insofar as they predict intrinsic ellipticity 
correlations of 1 to 10 % of the lensing signal for surveys with 
median redshift around unity, whereas the effect may even 
dominate the cosmic shear signal for shallow surveys. This 
behavior is expected because a high median redshift corresponds 
to a broad redshift distribution of galaxies, which decreases 
the probabil ity of finding two galaxies at the same redshift. 
iBrown et al I (|2002|) detected intrinsic ellipticity correlations in 
the SuperCosmos field, which is shallow, so that all measured 
ellipticity correlations h ave to be purely intrinsic; see also 
iMandelbaum et alj (l2006h for an analysis of SDSS data. 

So far, only toy models exist for the i ntrinsic alignment 
of galaxies (see e.g. iHirata & Seliafl |2004 HS04 hereafter), 
since little is known about the underlying processes of galaxy 
formation and evolution in their dark matter halos. However, 
as galaxies need to be physically close to mutually align, 
intrinsic ellipticity correlations are relatively easy to remove 
from the cosmic shear signal if red shift information is avail- 
able. iKing & Schneided (l2002i l2003h and lHevmans & HeavensI 
(I2OO3I) used different methods for downweighting galaxy pairs 
close in redshift, thereby reducing contamination by intrinsic 
alignment. Besides, intrinsic ellipticity c orrelations can gener- 
ate B-modes (HS04. lHevmans et al.ll2006h . which could be used 
for its identification if alternative sources of a curl-component 
can be excluded. 

The third term in ([T]l describes the correlation between the 
intrinsic ellipticity of a background object and the shear signal 
of a foreground galaxy. As the lensing effect is generated by the 
matter structure between z = and z = z,, there is no physical 
interaction between the matter acting as the lens and the galaxy 
at Zj, so that the correlation is expected to be identical to zero. 

The last term can yield a contribution if the matter struc- 
ture generates a tidal gravitational field in which a close-by 
galaxy aligns and by which a background object is gravitation- 
ally lensed. The foreground and background galaxies are found 
to be primarily oriented perpendicular to each other, resulting in 
an anti-correlation of the ellipticities. On small scales, a correla- 
tion might even ste m from the alignment o f a foreground galaxy 
with its own halo dBridle & Abdallall2007l) . We refer to this ef- 
fect, first mentioned by HS04 as 'shear-ellipticity correlation' 
(although a variety of names exists in the literature). 

In contrast to the foregoing case, however, the removal of 
shear-ellipticity correlations is more complicated, as they are 
not restricted to physically close objects. On the contrary, their 
contribution increases for larger separations in redshift, becom- 
ing the dominant contamination for deep surveys (HS04). N- 
body simulations yield an upper limit of about 10 % for shear- 
ellipticity correlations (Heymans et al. 2006), in agreement with 
the analytical approach of HS04. Such a contamination yields 
an underestimation (due to the anti-correlation) of erg of about 



5 % and a bias on the parameters of the dark energy equa- 
tion of state by up to 50 % if acting together with intrinsic 
e llipticit y correlations (Bri dle & King 2007, BK07 hereafter). 
iMandelbaum et alj (l2006l) report the first observational verifi- 
cation of shear-ellipticity correlations in the Sloan Digital Sky 
Survey (S DSS) and estimated an upper limit of 20 % contamina- 
tion, while iHirata et alj d2007D find a best-fit intrinsic alignment 
model with data from SDSS and the 2SLAQ survey, which pre- 
dicts 6.5 % contamination to the cosmic shear signal for a survey 
of similar galaxies. 

The methods suggested for eliminating shear-ellipticity cor- 
relations from the cosmic shear signal mostly take advantage of 
the different red shift depend ence compared to a pure lensing sig- 
nal (e.g. HS04. lKing|l2005L BK 07), but they sti l l partl y rely on 
the uncertain models. Although ^ Hevmans et al.l d200 6^ reported 
the detection of weak B-modes caused by shear-ellipticity cor- 
relations in their simulations, it is questionable whether this ef- 
fect is strong enough to be unambiguously measured and there- 
fore apt to remove the contamination. To provide a secure means 
of arriving at a data set free of shear-ellipticity correlations, we 
propose a purely geometrical method in the following, applica- 
ble to a cosmic shear survey with photometric redshift informa- 
tion. Making use of the characteristic dependence on redshift of 
these correlations, as already pointed out by HS04, we eliminate 
('null') the contributions to the cosmic shear signal that are pos- 
sibly subject to a contamination by intrinsic alignment. 

In Sect.m we introduce o ur notation, la rgely based on 
iBartelmann & SchneideJ d2001h and lSchneideJ (12006.) . and out- 
line the basic elements of our method. We present two numer- 
ical and one analytical approach to determining weight func- 
tions, modifying the lensing efficiency to achieve the nulling of 
signals, in Sects. [TTI to 13.31 Section [3741 provides details on the 
implementation and the results for these approaches, while in 
Sect. |3.5| the concept of weight functions is generalized to higher 
orders. The loss of information caused by the elimination of 
shear-ellipticity correlations and its consequences on constraints 
of cosmological parameters is discussed in Sect.H] Finally, we 
conclude in Sect.|5] 

2. Method 

The gravitational lensing of distant sources by the large-scale 
structure can be represented by the deflection of light in a sin- 
gle equivalent lens plane, neglecting lens-lens coupling and us- 
ing the Born approximation. Accordingly, the three-dimensional 
mass distribution is projected onto this plane along the line of 
sight. Hence, for every comoving distance x of a source, the di- 
mensionless surface mass density or convergence k{0,x) can be 
obtained from the three-dimensional density contrast 6 via 



K{e,x) 



Jo Mx) a(x') 



(2) 



2c2 Jo Mx) 
where a(x) is the cosmic scale factor and 

I sin( 

Mx) ^\x K^O (3) 

[ sinh( yTKx)! K<Q . 

For more details see e.g. [Schneider! d2006l) . In the following, the 
Universe is assumed for simplicity to be spatially flat, i.e. K -0. 
If photometric redshift information is available, several projec- 
tions of k{0,x) can be constructed. 



/^«(6») 



r 

Jo 



dxp"\x)<0,x) 



(4) 
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where the upper hmit of integration in (|4| is the comoving dis- 
tance horizon. The weight function denotes the normal- 
ized probabiHty distribution of distances for a galaxy population 
In this work we assume that these distributions are disjunct 
redshift bins, i.e. we neglect the errors of photometric redshifts. 

Consider now a tomographic shear signal from redshift bins 
/ and j, obtained via a second-order cosmic shear measure such 
as the coiTelation functions ^'^^\6), define d as usually in th e con- 
text of weak gravitational lensing (e.g. ISchneidenuOOq) . The 
contamination by shear-ellipticity correlations is caused by the 
matter distribution within the redshift slice situated closer to the 
observer, see ([T]) and the discussion thereafter. If one cannot rely 
on a model governing the statistical properties of this contami- 
nation, it is necessary to eliminate the contribution of the matter 
within the lower redshift bin to the cosmic shear signal com- 
pletely because only then can one ensure that shear-ellipticity 
correlations are not present anymore. Thus, the weighting in the 
projection of the density contrast has to be modified appropri- 
ately. Instead of using the galaxy distance distribution as in (|4|, 
we now write 

K^'\e):= dx B^'>(X) K(e,x) , (5) 

Jo 

where B^'^(x) denotes the new weight function, measured in units 
of inverse length like p^'^ix)- After inserting (|2]l into ^ and re- 
arranging the integrations, the relation of k{6) to the density con- 
trast reads 

\0) = - , Ax — — — 5 (xO,x) , (6) 

2c^ Jo a(x) 

with the lensing efficiency 

1*%) = jr''"dyB<'V)(i-^) . (7) 

The lensing efficiency determines the amplitude by which the 
density contrast at the coiTesponding distance contributes to the 
convergence. Thus, demanding that the matter structure at a co- 
moving distance does not yield any contribution to the cosmic 
shear signal is equivalent to the constraint 

f^(Xi) = jT^" dx (1 - = . (8) 

Equation (O ensures the 'nulling' of the contribution of mat- 
ter at Xi to the shear signal, and consequently of the shear- 
ellipticity correlations, motivating the name 'shear-ellipticity 
nulling' of this method. It makes use of the characteristic de- 
pendence on distance, 1 - Xi/x, which corresponds to the ra- 
tio Dis/D^ of the angular diameter distance between lens and 
source and the one between observer and source. The nulling 
te chnique has been a pplied before in the context of cosmic shear 
bv lHuterer & Whitd (I2005h who aimed at nulling small-scale in- 
formation influenced by poorly determined baryonic physics. As 
Dds/Ds is a smooth function of both the lens and the source po- 
sition, the integration in (|8j is rather insensitive to small changes 
in the distances of the mass distribution, acting as the lens. 
Consequently, the contributions from distances slightly smaller 
or larger than;f^, are strongly downweighted by B^'^ix), too, pro- 
vided that the weight function is also smooth. This leads to the 
near cancellation of the signal from the whole bin / ifxi is chosen 
to be the distance corresponding to the center of the bin. 

For reasons of simplicity, we assume in the following that 
the binning in redshift is equidistant with width Az. However, 



the subsequent equations can also be formulated using variable 
bin sizes. We define a new second-order measure in terms of the 
tomography correlation functions as 

H«(0):= dxB'^'\x)^Az(Xi),z(x),e), (9) 
Jo 

where the correlation functions ^±(zi ,Z2,ff) are - formally - eval- 
uated at the exact, e.g. spectroscopically determined, redshifts zi 
and Z2- So far, B^'^(x) is only constrained for xt < x ^ Xhoi by 
(|8j. We set B^'^(x) = for < x ^ Xi because otherwise cor- 
relation functions with z(x) < zixd, which can be contaminated 
by intrinsic alignments at z(x), would contribute to (|9]l. Through 
this, no information is discarded since the correlation functions 
are symmetric in their first and second arguments, so that, swap- 
ping the redshifts, they enter other measures, constructed like in 
(|9j but with smaller ;f^,. 

Now assume that ^±(z(Xi), z(x), S) is contaminated by shear- 
ellipticity coiTelations. Since B^'^{x) vanishes for x < Xi^ it is 
sufficient to consider z(x) > z(Xi) in the integral. Then intrinsic 
alignment has to be generated by the matter structure at;f', ; how- 
ever, the contribution of matter atXi to the convergence and sub- 
sequently to the shear entering the correlation function is elim- 
inated by the weight function when chosen according to dHJ, so 
that B}+(6) is free of shear-ellipticity correlations. 

Since in practice only photometric redshift data will be avail- 
able, we first transform (|9]l to z as the integration variable and 
subsequently approximate the integral as a Riemannian sum, 

H«(0) ^ B^'^Wzj)) &\0) X'izj) Az, (10) 

.1=1 

where A^- is the number of redshift bins andx'(z) the derivative 
of the comoving distance with respect to z, calculated from 

Xiz) = r dz' [njl + z'? + ^aV"^ . (11) 
m Jo 

The distance-redshift relation enters into many of the following 
equations since (|2|l and related equations are formulated in terms 
of comoving distance, whereas the binning is performed in red- 
shift, which is the observable. Henceforth, we refer to the first 
bin of the tomography correlation functions entering (fTOl i. i.e. 
the bin where the signal is nulled, as the 'initial bin'. In ( flOl ) 
Zj denotes the lower boundary redshift of bin j, implying that 

in the sum the correlation function ^^i^\6) is registered onto the 
lower boundary of bin j. To render ( fTOl l a good approximation, 
the redshift binning has to be sufficiently narrow. In addition, it is 
assumed that the redshift bins cover the total galaxy population 
almost completely, so that on discretizing (|9j, the upper bound- 
ary of the integral can be reduced to the upper boundary of the 
highest redshift bin, denoted by z^ax- 

For the sake of lower computational efforts, the further in- 
vestigation of shear-ellipticity nulling is performed in Fourier 
space, by defining the new power spectrum 

rXhi,, 

n«(f) := dxB^'>(x)PAz(Xi),z(x)J) 
Jo 

^ Y^B^\xizj))PV\()x'izj)Az, (12) 

where Pk{z\,Z2, () denotes the convergence power spectrum for 
exactly known redshifts, and the angular frequency £ is the 
Fourier variable on the sky. The quantity P^i^''\{) constitutes 
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the corresponding tomographic measure, related to the three- 
dimensional power spectrum of matter density fluctuations Pg 
via Limber's equation in Fourier space (Kaiser 1992), 



4c4 



Jo 



(13) 



where the lensing efficiency g^'^Cy) in its original form, i.e. with 
p'''^(x) as weight, enters. However, since shear-ellipticity nulling 
requires detailed redshift information, Az <K 1. The weighted 
convergence in (HI can then be written as k^'^(B) k, K{6,Xi), or 
equivalently as p^'^ix) ~ SD(x~Xi) to good approximation, where 
6u is the Dirac delta distribution. In our case the dependence of 
the tomography power spectra on the probability distribution of 
galaxy distances is only marginal. 

From the N,(N, + l)/2 tomography power spectra, which are 
in principle available, only A^. power spectra n*''(^') are con- 
structed via (fT2l i. However, we determine more than just one 
weight function B^'^(x) per initial bin. Thus, as explained in 
Sects. [33] and|4l the number of new power spectra will be equal 
to the number of input P^i^\{). Note that, after the explicit com- 
putation of weight functions, we show newly constructed power 
spectra and discuss their features; see Sect. l3.5l 

As the weighting by B^'^(x) does not depend on angular 
scales, the considerations made above for also hold for 

n'^''(^). The weight functions, once obtained, can be applied 
directly to measures accessible to observations, constructed in 
analogy to ( fTOl i. by means of the transformations between con- 
vergenc e power spectra and rea l-space measures, as for instance 
given in [Schneider et alj (l2002h . Hence, ^ and ([T2] ) are related 



via 



S®(0) 



Jo 27r 



(14) 



where denotes the Bessel function of the first kind of order /i. 

In ([l2b it is sufficient to let the sum over redshift bins start 
only at ! H- 2 because all other function values of B^'^(x) vanish. 
First, we demanded that B*''Ct) = for Q < x ^ Xi-^^ addition, 
to avoid contamination by intrinsic ellipticity correlations, we 
force the weight function to vanish within bin / completely, im- 
plying that B'^^^ix) departs from zero only from the upper bound- 
ary of bin /, i.e. from redshift Zi+\, onwards, so that the first non- 
vanishing value of the weight function in ([T2l l is Z?*'*(^(z,+2)). 

Consequently, P^i^\{) with \i - i\ < 1, i.e. power spec- 
tra which auto-correlate redshift bins or cross-correlate adjacent 
bins are removed from the new measure ( [T2] l. se curing the down- 
weigh ting of intrinsic ellipticity correlations. iTakada & White! 
( |2004|) demonstrate that errors on inferred parameters degrade 
by only about 10% for 5 or more redshift bins if the auto- 
correlation power spectra are simply excluded from the analy- 
sis in shear tomography, which constitutes a simple and efficient 
method. Due to the narrow redshift binning our method requires, 
the contamination by intrinsic alignment may extend to neigh- 
boring bins (see e.g. BK07). Since we assume a choice of the 
redshift bins such that the number of galaxies with z > Zmax 
is negligible, we set B^'^(x) = for x ^ xi^max)- The conven- 
tions concerning the construction of weight functions made in 
this section are summarized in Fig.[T| 

3. Determination of weight functions 

To show the feasibility of shear-ellipticity nulling, it is impor- 
tant to know to what extent and how efficient the weight func- 
tion B^'^ix) can be constructed in practice. As a basic condition. 



(ij=i+2) (ij=i+3) 
P„ P, 




<ij=N,-l) (ii=N,) 
P P 



N,-l N\ bin no. 



Fig. 1. Sketch illustrating the conventions made for the con- 
struction of weight functions. The redshift Zj denotes the lower 

boundary of bin j; the tomography power spectrum P^!,^\() 
(or the corresponding real space measure) is registered onto 
this redshift in ( fT2] i. With respect to these boundary redshifts, 
B^'HxiZi+i)) is the first and B'^'HxizN,)) the last non-zero value of 
the weight function. The upper boundary redshift of the highest 
bin is denoted by Zmax- The initial bin is located between Zi and 
Zi+\ with a central redshift of zi, so that^i - xizd- 

B^'Hx) has to obey the constraint ([8]l. In addition, it is required 
that the weight function should be optimized in the sense that, 
when aiming at constraining cosmological parameters, the in- 
formation content in the new power spectrum W^'HE) attains a 
maximum. This condition is quantified in terms of the Fisher 
matrix, which reads, when measuring the power spectrum in N( 
angular frequency bins, as 



f(0 



Z 

a, 13=1 



5n«(4) 



f \ 'o-/ 



dpv 



(15) 



The vector p is composed of the cosmological parameters that 
are considered. Its size A^^, implies the dimension A^^, x Np for 
the Fisher matrix. The index f is assigned to the power spec- 
trum n^'*(^), indicating that it is evaluated at the parameters of 
the fiducial cosmological model. The covariance C^'' measures 
the correlation of n*^'' between different angular frequency bins. 
From the definition (fTST i one finds 



k=i+2, /= j+2 

X B^'^ixizk)) B^jHx(zi))x'(zk)x'(zi) Az^ , 



(16) 



where the covariance of the convergence power spectra is given 
by (seelJoachimi et al. 2008 for details) 



{cr\ 



In 



witii P 



m ._ p(ki) 



pr + sm 



2n« 



(17) 



where A denotes the size of the survey and /S.i the width of 
the angular frequency bins. The shape noise contribution is 
governed by the intrinsic galaxy ellipticity dispersion cr^, i.e. 
crl - (eV") and by the mean number density of galaxies 
in bin /. Note that for ( [TtI i to be valid, we have assumed that the 
shear field is Gaussian and that the survey geometry is simple. 
Consequently, the power spectrum covariance is diagonal in an- 
gular frequency space, so that ( fTSI l is diagonal, too, rendering the 
inversion needed for ( fTST i trivial. 

Besides, in ( [TtT i it is assumed that any B-mode contribu- 
tion to the power spectrum vanishes, which is expected for a 
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pure lensing signal and which was already implicitly assumed 
throughout the preceding derivations. As a side remark, if there 
was a B-mode contribution due to shear-ellipticity correlations, 
the nulling technique could eliminate it as long as the B-mode 
signal has the same characteristic redshift dep endence as the 
E-mo des. However, though seen in simulations (iHevmans et alj 
I2OO6I) . B-mode shear-ellipticity correlations are not expected 
from theory (HS04), so that the scaling with redshift is unknown. 

As the actual quantity that is maximized, we choose the trace 
of the Fisher matrix because of its simple functional form, with 
the Fisher matrix elements entering tr(F^'*) linearly. Actually, 
the determinant of the Fisher matrix would be a more natural 
choice since its inverse is a measure of the volume of the error 
ellipsoid in parameter space, the inverse Fisher matrix being an 
estimate of the corresponding covariance matrix. However, due 
to the non-linearity of the determinant the numerical treatment 
is considerably less stable. As we demonstrate in Sect.lH the 
trace fulfills the demand of concentrating the bulk of information 
into the power spectrum constructed out of the optimized weight 
function well enough. Exemplary, lower-dimensional calcula- 
tions show that the difference in the form of the weight functions, 
using the trace or the determinant, is marginal. 

Considering the continuous limit of Tl^'\{) in ( fTSl i for a mo- 
ment, one arrives at a problem of variational calculus, determin- 
ing the maximum of tr(F*'') with respect to the function B^'^(x) 
under the constraint However, due to the binned redshift in- 
formation, B^'^ix) only enters the equations in the form of dis- 
crete function values, turning the problem into a maximization 
with respect to these function values. The complete weight func- 
tion is then constructed by either interpolation or in parameter- 
ized form. Plugging ( fT6b and ( fTTI ) into ( fTSl l yields 



N[ 



(18) 



j,k=i+2 




■XiZj)+Pri(a)- 



B(x(Zj)) 



All power spectra and derivatives thereof are evaluated at the 
fiducial values of the cosmological model. Note that tr(F''^) is 
independent of the overall amplitude of as must be the 

case. The function values of B^'^(x) enter tr(F*'^) non-linearly, so 
that analytical progress is hampered in the general case. Both 
numerical approaches and analytical approximations are investi- 
gated in the following. 



3. 1 . Piecewise linear approach 

First, we consider a piecewise linear ansatz and write 

B(%(z)) Bj + {Bj^i - B,] for z e [zj, Zj^i] , (19) 

where the notation Bj = B^'\xizj)) was introduced for conve- 
nience. Moreover, we identify B^^+i = B^'\xiZmdii)) - 0. The 
superscript (/) of the weight function is dropped when using this 
shorthand notation. As before, zj denotes the lower boundary 
redshift of bin j, whereas we denote the central redshift of the 



initial bin, entering the constraint (jSj, as z,-. With the weight func- 
tion in the form of (fT9] l, the constraint reads 



j=i+i 



y(^)(l_^J = 0.(20) 



Due to the conditions imposed on B^'^ix), as mentioned in 
Sect.|2] - holds, so that the constraint assumes the com- 
pact form 



N- 



j=i+2 

where the quantities 



(21) 



(22) 



= X(Zj+i) -x(Zj) -Xizi) In 



n*' z-z, , 

J dz x (z) 



Az 



x(zj+\) 
X(Zj) 
Xizj) 
X(z) 



were introduced. Fixing one of the function values via ( 1211 1. the 
remaining Bj are varied to obtain a maximum of the trace of the 
Fisher matrix, using a Nelder-Mead simplex algorithm, which 
is robust and does not require partial derivatives. However, this 
type of maximization routine easily gets stuck in local extrema, 
which in this case are caused by single outliers among the Bj. 
These are considered unrealistic since B^'^(x) is expected to be 
smooth, but the piecewise linear ansatz does not put any con- 
straints on the derivatives of the weight function. 

To avoid the false maxima due to outliers, we subtract the 
term 



U:-AY,{Bj.i-2Bj + Bj^,) 



2n 



(23) 



>='■+! 



from ( fTSl l, summing the differences in slopes at the nodes, 
thereby disfavoring solutions with abrupt changes in the first 
derivative. The fudge factor A > can be chosen arbitrarily, 
while denotes a small, positive integer. The larger 77, the more 
weight is given to large changes in slope in U, where we found a 
suitable value of 77 = 2. If a value of A is chosen such that U and 
ti{F^'^) are roughly the same order of magnitude, the resulting 
weight functions are very smooth. By gradually lowering A, less 
smooth B^'^ix) with a more pronounced maximum are obtained. 
In case an outlier Bj occurs, the initial values of the simplex al- 
gorithm are altered, until a stable solution with A = results. 

The weight functions of this section are hardly susceptible 
to false maxima due to their smoothness, so that mostly we can 
set A = from the beginning. However, the procedure outlined 
here is necessary for the more oscillatory weight functions that 
will be computed in Sect. |3.5l Still, the final results presented in 
this work have been obtained with A = throughout. As can be 
seen from ( fTSl l, the maximum of tr{F^'^) does not depend on the 
overall amplitude of B^'^{x), leading to a degeneracy in the max- 
imized Bj, which will be hfted by a normalization, see Sect. l3.4l 

3.2. Chebyshev series approacii 

The second numerical approach assumes that B<'%) IS com- 
posed of a finite series of ansatz functions with a set of free 
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parameters. In this case, we choose Chebyshev polynomials of 
the first kind T^,, which already lead to good approximations for 
low polynomial orders and yield evenly distributed errors. The 
weight function is expanded as 



(24) 



f=0 



2y- kfc+i) +x(zmiLx)} 



where the b/j denote the A^c + 1 free coefficients. The argument of 
is chosen such that it takes on values in the interval [-1,1]. 
Plugging this definition into (O, one gets 



No 

with the definition 



dz Ixiz) - X(Zi+l )) Ixiz) - X(Zmax)} 



(25) 



(26) 



XiZm^!i) - X(Zi+\) 



Xiz) 



Again, one of the parameters is fixed by dZST l, while Nc pa- 
rameters are used for the multi-dimensional maximization of the 
trace of the Fisher matrix with the simplex algorithm. If the num- 
ber of free parameters is chosen to be more than about 5, the 
resulting weight functions are prone to significant oscillations, 
generated by the response of the ansatz polynomials to the steep 
rise in B^'^(x) foTx(z) - xizi+i) ^ 1- These unphysical features, 
corresponding to shallow maxima in parameter space, are readily 
detected by visual inspection of the resulting weight functions 
and avoided by altering the - in this approach low-dimensional 
- set of initial values fe^. As for the piecewise linear ansatz, we 
observe the degeneracy in the parameters yielding the maximum 
due to the free scaling of the weight function. 

3.3. Simplified analytical approach 

We elaborate on an analytical approach that is computationally 
fast and can provide an important consistency check for the pre- 
ceding numerical methods. However, as stated above, the non- 
linearity hinders the analytical treatment of the full problem; 
instead, we confine ourselves in the following to considering a 
single angular frequency bin and only one element of the Fisher 
matrix, i.e. a single cosmological parameter A vector notation is 
introduced as follows. 

Let the non-vanishing values of the weight function B j, given 
in their shorthand notation of Sect. 13. II form a vector B. Note 
that in the vector notation we again drop the superscript (/) since 
the initial bin that B refers to will be clear from the context. By 
defining another vector / with components 



X(zi) 
Xizj) 



X'izj) for = / + 2, 



(27) 



the constraint (O simply turns in its discretized version into 
(B-f)^O. (28) 

Note that the constant bin width Az is not included in /, as it 
drops out when setting the discrete constraint expression to zero. 
The covariance of n''\ given by (fTSI l. reduces to a scalar quantity 



due to the single angular frequency bin ( under consideration. 
Defining a matrix C with elements 



Cu:-{cf"\ x'(zk)x'iz,), 
the covariance can be written as 



(29) 



(30) 



C[;" = B^CB Az^ . 

The Fisher matrix element, now indicated by a subscript o, reads 
as 

i^® = (31) 
B^CB 

with a further vector defined for convenience. 



PJ ■-- 



(32) 



dp dp 

By means of this vector, one is able to rewrite the derivative of 
the power spectrum with respect to the remaining cosmological 
parameter p as dll^'\{)/dp - (B ■ p) Az, which can be seen by 
taking the derivative of ( fT2b . The bin widths Az cancel in Fg\ 
so that they do not need to appear in the definitions of C and p. 
The constraint is incorporated by means of a Lagrange multiplier 
A, leading to the quantity G :- f'"^ + A(B ■ /), which is to be 
maximized with respect to the components of the vector B. One 
obtains 



-2CB 



B^CB \B^CB, 
which can be formally solved for B, resulting in 



B = NC 



B^CB 



{B-p)^ 2 ^ \{B-p) 



B^CB 



(33) 



(34) 



where the free normalization N of B has been introduced. 

An overall scaling factor in B neither modifies the informa- 
tion content of 11*'* nor does it alter the constraint (O, which 
illustrates that the conditions stated above do not fix the normal- 
ization. The formal solution is plugged into ( |28] ), which is then 
solved for the Lagrange multiplier. 



{B-p) fC-'p 



2 B^CB fC-^f 
Replacing Aj2 in (l34l l subsequently yields 
B^CB 



B ^ N 



(B-p) 
N'ic'p 



fc-'f 



(35) 



(36) 



fc-'f 



where the scalar quantity {B ■ p) /B^CB was absorbed into the 
normalization N', so that now the righthand side depends no 
longer on B. Optimized analytical weight functions can be cal- 
culated by means of ( |36] |. interpolating linearly between the 
values of the components of B. In order to achieve results as 
close to the non-simplified, numerical approaches as possible, 
the employed element F*'* is chosen to be the diagonal element 
of the Fisher matrix that yields the largest contribution to the 
trace. Afterwards B is determined on a grid of ^'-values within 
the range considered in the numerical approaches, the solution 
vector resulting in the largest F*,'' being taken as the 'optimal' 
weight function. As long as this term, which is supposed to be 
the strongest contribution to the sum in ( fTSl l, dominates the trace 
of the Fisher matrix, we expect this procedure to yield a reason- 
ably good approximation to the numerical results. 
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Fig. 2. Weight functions as a function of redshift for A^^ = 40. The initial bins are given in the respective panels. Plotted are the 
simplified analytical solution as solid curve, the Chebyshev series solution as dashed line, and the piecewise linear solution as dotted 
curve. 



3.4. Setup & results 

To construct actual weight functions B^'^(x), a Active tomo- 
graphic cosmic shear survey with a comparatively large number 
of narrow redshift bins is needed. Due to the choice of n*'*(^) 
as the quantity considered, the necessary input data comprises 
a set of tomographic power spectra, which are obtained for a 
ACDM universe with fiducial parameters Qm = 0.3, Qa - 0.7, 
and Hq = 100/;i()okm/s/Mpc with /iioo - 0.7. The three- 
dimensional power spectrum of density fluctuations is specified 
by the primordial slope - I, the normalizat ion erg - 0.9 and 
the shape parameter F, calculated according to Sugivamal ( Il995h 
with Qb - 0.04. The line ar power spectrum is given by the fit 
formula of iBardeen et al.l 0984), w hile the non - Unear evolution 
is included via the prescription of ISmith et alj (l2003h . The to- 
mography power spectra are then determined for N( - 75 loga- 
rithmic angular frequency bins between ^ = 50 and i = 10**. 

Furthermore, we specify survey properties that enter ( fTSl l 
via the power spectrum covariance ( fTTl ). For this a normalized 
galaxy redshift probability distribution 

with zo = 10 and /3 = 1.5 is assumed jSmail et al.lll994l) . 
The redshift distribution is cut off at z^ax = 4, requiring a 
renormalization, which leads to the modified distribution Pcutiz). 
However, due to the large cut-off redshift the modification is 
marginal. The dependence of the covariance on the survey size 



A is trivial, its value being irrelevant for the determination of the 
weight functions. For later calculations of likelihoods, we set A 
to a fiducial size of 1 deg^. Moreover, we set the intrinsic ellip- 
ticity dispersion to <t^ - 0.4 and choose a mean galaxy number 
density of n = BOarcmin"^. The bin-wise number densities are 
obtained by 

«('■> = h dz /,eut(z) . (38) 

The derivatives of the power spectra with respect to cosmologi- 
cal parameters in ( fTST l are obtained via finite differencing, while 
dx'{Zj)ldpf, is calculated from (fTTT i in analytical form. 

To determine the Fisher information, we considered the set 
of cosmological parameters p - (Qm, Qa, (Xg, /iioo). Aiming at 
smooth weight functions, we used a large number of redshift 
bins, i.e. A^, = 40, corresponding to Az - 0.1. As mentioned 
in the foregoing section, the normalization of the weight func- 
tions is not yet fixed. To allow for direct comparison of the three 
approaches, we imposed the condition 

Ax |B®Or)l' = 1 . (39) 

•Jxi 

Furthermore, the free sign of B^'\x) is chosen such that the 
weight function first assumes positive values when departing 
from zero at the upper boundary of the initial bin. 

In Fig.|2] the resulting weight functions of all three ap- 
proaches under consideration are shown for varying initial bin 
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Fig. 3. Weight functions as a function of redshift for two differ- 
ent redshift binnings. Solid curves were obtained with A^. = 40; 
dashed lines with = 20. Plotted are the simplified analytical 
solutions for initial bins 2, 6, and 10 in the case of 40 bins, and 
1, 3, and 5 in the case of 20 bins. 



i. All methods are in very good agreement; only in the upper 
left panel are larger deviations visible. A close inspection reveals 
that the weight functions constructed by means of the Chebyshev 
series suffer from slight oscillations that can rapidly increase in 
amplitude in some cases if A^^, is chosen too large. These can pre- 
sumably be explained by the steep rise of B^'^'ix) near the initial 
bin. Generally speaking, the agreement justifies the assumptions 
made in the different approaches; in particular, the results of the 
relatively crude approximations of the analytical approach are 
compatible with the numerical calculations, so that it is well- 
suited to further investigation. The weight functions have a zero- 
crossing, which is expected due to where the term 1 - Xilx 
is non-negative throughout the integration interval. The largest 
weight is assigned to those redshifts that are located directly 
above the initial bin, respectively, because the efficiency of the 
lensing of a source in this range of redshifts by the mass distri- 
bution within the initial bin is low or, in other words, Dds/£>s is 
small, which decreases the contribution by shear-ellipticity cor- 
relations. 

The division of a cosmic shear survey into 40 redshift bins 
is realistic in the near future; however, the bin size would not be 
chosen constant as in this study, but probably scale with 1 + z. 
Mainly for computational reasons, we reduce the number of red- 
shift bins used in the following likelihood analysis to 20, a num- 
ber achieved by some of the upcoming wide-field projects such 
as Pan-Starrs, KIDS/Viking, the Dark Energy Survey, LSST, or 
Euclid. We compare the form of the weight functions obtained 
above with an analogous set, determined for A^, = 20. Figure [3] 
illustrates for a sample of analytical solutions for B^'^(x) that the 
sets for both redshift binnings agree well. The higher density of 
sampling points in the case of A^. = 40 enables a steeper rise 
of the weight functions at the upper boundary of the initial bin, 
leading to deviations in similar magnitude in the tail of B^'^(x). 
These results also suggest that the effects due to the discretiza- 
tion of the weight functions are negligible as long as the number 
of redshift bins is not chosen too small. 



3.5. Higher orders 

The power spectrum Yl^'Ht), given by (fT2l i. is a linear combina- 
tion of the convergence power spectra, with the weighting de- 
termined by the function as calculated in the foregoing 



section. More such linear combinations can be constructed with 
differing weight functions that still obey the constraint equation, 
resulting in further power spectra free of shear-ellipticity corre- 
lations. If one retains the condition of maximizing the Fisher ma- 
trix and, in addition, demands that the weight functions should 
be orthogonal with respect to each other in a suitably defined 
sense, one arrives at higher-order measures that have the second- 
most, third-most, etc., information content. 

Since in ( fT2] i the first term that yields a contribution is for 
i - i + 2 (see also Fig.[T]), N~-i-l convergence power spectra 
are used to form in the implementation presented above. 

As a consequence, one is able to construct N,-i-2 mutually or- 
thogonal power spectra Tl^'\£) from this data set. The additional 
combination that could be built furthermore with linear indepen- 
dence from the set of convergence power spectra then necessar- 
ily violates ([8]l and consequently contains shear-ellipticity corre- 
lations. 

Denoting the order of the weight function by a subscript in 
square brackets, the condition of mutual orthogonality between 
weight functions of order q and r can be formulated as 



•A-(Zn,ax) 



(40) 



for all orders q > r, where w{z) is an arbitrary weight func- 
tion. As far as the two numerical approaches are concerned, the 
higher order weight functions are obtained by fixing one further 
free parameter for every orthogonality condition. In the case of 
the piecewise linear ansatz, one obtains by plugging the ansatz 
functions (O into ( l40l l 



,/=l+l ^'i 



dz 



Z-Zj 



[qlj+\ - B[q]J 

(0 



X 



B' '(xiz)) w(z) = (41) 



for every r = 1 , .. , <7 - 1 . The lower order weight functions have 
been determined in advance and are known. Defining 



dz B{;.'Or(z)) w(z) , 



(42) 



^ - r 

one can write (HTt in analogy to (|2TI) as 



dz ^ Blpcizy) w(z) , 



The Chebyshev approach yields, inserting (|24] | into (|40] |. 



(43) 



(44) 



again for r - 1, .. , ^ - 1, where we have defined 

:= V^'^^ Az\x{z)-xizM)\\x{z)-x{z^..)\ (45) 



/ 2^(^) - \xiZi+\) +;r(Zmax)) 



I ^'(Zniax) -^(Z/+l) / 



B^Mz)) w(z) . 



Besides, higher order weight functions can be constructed 
with the analytical ansatz, again considering only a single angu- 
lar frequency bin and solely one component of the Fisher matrix. 



B. Joachimi and P. Schneider: The removal of shear-elhpticity correlations from the cosmic shear signal via nulling techniques 



9 





1.5 




1.0 




0.5 








0.0 




-0.5 




1.5 




1.0 




0.5 


X 






0.0 




-0.5 




1.5 




1.0 




0.5 








0.0 




-0.5 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

initial bin: 2 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I i_ 
initial bin: 4 




Fig. 4. Second-order weight functions as a function of redshift for N, - 40. The initial bins are given in the respective panels. The 
coding of the curves is the same as in Fig.|2l 



optimizing it subsequently as outlined in Sect. 13. 31 To clarify the 
notation, we rewrite the result for the first-order weight function 
(l36b as 



B 



[1] 



with p 



fc-'f 



f 



(46) 



As in Sect. 13. 31 we perform the derivation for an initial bin /, 
where the index does not explicitly appear in the formulae, but 
enters the quantities /, C, and p. In this context the condition of 
orthogonality can be implemented as 



(47) 



for all orders q > r, where Bj := B'^'^(x{zj))w{zj) was defined. 
This condition is incorporated into the maximization by more 
Lagrange multipliers, the expression to be maximized for order 
q turning into 



9-1 



(48) 



where the yu^^^j are the Lagrange multipliers for the respective 

orthogonality conditions. The Fisher matrix element F*'' is still 
given by ( |3TI ). now with weight functions B[q]- After taking the 
gradient with respect to the components of B[g] in analogy to the 
first-order calculation, one arrives at the formal solution 



= NC 



2 •' 



\2 



(49) 



The Lagrange multipliers are successively replaced by inserting 
this solution into the corresponding constraint equations. As a 
first step, from (j' ■ B[q]j - one obtains 



A 



[9] 



(50) 



Plugging in this expression, ( |49] | turns into 
9-1/;!-. ,B],CBu 



B 



[9] 



[ 7~i ^ ■P) 



with the definition 
fC'f 



f, 



(51) 



(52) 



where x e |p, B[cy]}. Again, multiplicative scalars have been ab- 
sorbed into the normalization. In a similar manner, inserting dSTT i 
into (l47T i for r - I leads to 



I r=2 ^ [Blc,]-P) 



(53) 
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where we set 



(54) 



If one continues likewise for the remaining constraint equations, 
one obtains in accordance with ( |46] | the compact result 



(55) 



for all orders q, where the redefined normalization is denoted by 
just N again. Here we made use of the recursion relation 



~{r-\] ^['■-1] 



(56) 



supplemented by the initial step ( |52] l. Consequently, a recursion 
relation for the vectors B^q^, corresponding to the higher order 
weight functions, can be derived, which reads 



B 



[«] 



N{B 



(57) 



for q > 2, the vector B[i] as the starting point for this recursion 
being given by ( |46] l. 

The weight function w{z) could for instance be chosen, such 
that it scales with the redshift probability distribution Pcutiz), 
assigning a larger weight to well-sampled redshift ranges. 
However, we set w(z) = 1 in the following for reasons of sim- 
plicity. In addition, (l47T l then turns into an orthogonality relation 
also for the vectors B[g]. 

In Fig.|4]the results for second-order weight functions fip^jCf) 
are plotted, for all three methods considered in this work and 
using the same setup as described in Sect. l3.4l Apart from slight 
numerical instabilities, as can be seen for the linear approach in 
the upper left panel, and the differing response to sharp peaks in 
the weight functions, most prominent in the lower right panel, 
the curves largely agree. Again, the simplified analytical ansatz 
proves to be compatible, being computationally advantageous 
to a large extent due to its recursive form ( ISTT i. For all curves 
the number of zeros has increased by two compared to the first- 
order results, a trend that continues for higher orders. A sample 
of analytical solutions for third- to fifth-order weight functions 
is given in Fig.|5] 

Adopting the notation introduced for the weight functions, 
the new power spectra can be generalized to higher orders as 



n« (^) : 



(58) 



for initial bin / and order q. Employing analytically determined 
weight functions, we compute these power spectra for different 
initial bins and orders, the resulting graphs shown in Fig.|6] For 
reference the convergence power spectrum Pk{{), integrated over 
the full redshift distribution as given by (IJTT i. is plotted in addi- 
tion. It is important to note that n|^^j(^ and P^it) can only be 
compared with difficulty in terms of the overall amplitude, since 
for the newly constructed power spectra, the amplitude can be 
chosen arbitrarily due to the free normalization of the B^^^-^(x). 
In Fig.|6] it is fixed by ( |39] l, so that the weights are of order 



unity. Therefore it is evident that the n|^*j(^) have considerably 



1.0 
^ 0.5 

-0.5 
1.0 

^ 0.5 



-0.5 h 
1.0 



0.5 - 



^ 0.0 
-0.5 



I ) 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

,* 3rd order initial bin 

''M — 2 : 

4 - 
6 
8 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



4th order 




■■■■■-7 
\ t 

1 1 1 1 1 I f If 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 hh' 



5th order 




-I I I I V i: I I I I I I 'i' I I \\ 







1 



2 
z 



Fig. 5. Higher-order weight functions as a function of redshift for 
- 40. From top to bottom the third to fifth order analytical 
solutions are shown. Solid curves correspond to initial bin no. 2, 
dotted curves to initial bin no. 4, short-dashed curves to initial 
bin no. 6, and long-dashed curves to initial bin no. 8. 



lower amplitude than the reference power spectrum since to ob- 
tain the former quantities, power spectrum signals are partially 
subtracted. 

Concerning shape, the n|^^j(/') show a largely similar behav- 
ior with respect to the convergence power spectrum, the latter 
peaking at higher values of i. This can be understood by tak- 
ing into account that the tomography power spectrum with the 
smallest difference between bins / and j contributes most to the 
respective n|^^j(^, as can be concluded from the form of the 
first-order weight functions, its pronounced peak being located 
just above the initial bin, see Fig.|2] Thus, the new power spec- 
tra receive their signal preferentially from less distant galaxies, 
so that they probe smaller physical separations for fixed angu- 
lar scale or I, respectively. On small scales non-linear structure 
evolution sets in, enhancing the signal. Consequently, the char- 
acteristic bump caused by non-linearity is visible for smaller 
i.e. larger angles, in the new power spectra in comparison with 
P^it), leading also to the shift of the peak. 
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Fig. 6. New power spectra as a function of angular frequency, 
making use of analytically determined weight functions. The 
power spectra n|^^j(f) are given as black curves, their order q 
ranging from 1 in the top panel to 3 in the bottom panel. Within 
each panel the power spectra for initial bins / = 1, ..,5 are 
plotted in the following sequence of line types: solid, dotted, 
dashed, chain-dashed, and dot-dashed. In addition, the conver- 
gence power spectrum P^iO, integrated over the redshift distri- 
bution (IJTI i. is shown for reference as gray curve. 



4. Information loss 

By eliminating contributions to the cosmic shear signal at cer- 
tain distances from the observer, one necessarily reduces the in- 
formation content of the data set, so that the desired constraints 
on cosmological parameters are less stringent. Thus, to judge 
the practical value of the nulling technique, we are going to 
quantify the accuracy with which cosmological parameters can 
be determined by the newly constructed power spectra (fTZb . If 
one considers the set of A^^ - ; - 1 tomography power spectra 
used to construct n|^^j(^) as the components of a data vector, 
then shear-ellipticity nulling is equivalent to a rotation of this 
vector such that all but one component of the resulting vector 



are free of shear-ellipticity correlations. The 'cleaned' compo- 
nents correspond to the N^ - i - 2 new power spectra (|58] l for 
q - 1 , .. ,Nz - i - 2, whereas the last component must contain a 
weight function that is coUinear to /, or simply 



N~ 

j=i 



('■;■) 



W 1- 



X(zj) 



X'izj) Az ; 



(59) 



i.e., / itself is chosen as the weight, see dZTl i. The vector rotation 
mentioned above is invertible, so that, illustratively, it is obvious 
that using the new, full data vector instead of the one containing 
the convergence power spectra for the data analysis, one should 
obtain the same results. This statement is equivalent to the Fisher 
matrix, as a measure of the information content, being invariant 
under such orthogon al transformations of the data vector (see 
lTegmarketalJll997h . 

Let the complete data vectors forming the basis of this anal- 
ysis be D for the original set and D' for the transformed one, 
which can be written in the convenient form 



D 



D' = (n|;j(A), ..,n|;^_2](^i)'n|ii(A), ..,ng_3j(A), 



n|^J(A), 



,n|^-^VA),n|^-^'(A),n|^-^'(A), 



n|;j(^2), 



>n[J'(^Aif), 



[1] 



(60) 



We refer to D as containing the full information although the 
vector is not composed of all tomography power spectra; how- 
ever, as already discussed above, these entries would most prob- 
ably have to be discarded anyway to avoid intrinsic ellipticity 
coiTelations. With the choice ( l60l l. both vectors have the same 
dimension A^d = N((N^ - 1){N^ - 2)/2. Their components are 
ordered such that the corresponding covariance matrices obtain 
a block-diagonal structure because the power spectra evaluated 
at different angular frequencies are not correlated due to the as- 
sumption of Gaussianity. Hence, the covariance Co of the data 
vector D reads 



' {Adi Ad]) 

{Ad2 Adl) 



(61) 



(Arf^, Arf]v,)j 

where for the sake of a compact notation, the vector 



,Pr^'(id,P?''({d, 



2N-) 



(A)) (62) 



was introduced, so that D = (d], .. ,dM^). The remaining non- 
trivial blocks {Adi A</[) for the angular frequency bin (i with di- 
mension (N^-l)(N^-2.)/2x(N,-l)(N:-2)/2 each are computed 
by means of ( fTTI l and can then readily be inverted numerically. 
The covariance of D' is dealt with analogously. 

Since the set of n|^'^(^) with q - I, .. ,Nr - i - I contains the 
full information, the optimization of the weight functions with 
respect to the trace of the Fisher matrix becomes superfluous in 
this situation. Instead, one can construct the B^''^ simply as a set 
of orthogonal vectors, starting with /, for instance by means of 
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Fig. 7. Contours of posterior likelihood before and after shear-ellipticity nulling. Shown are all possible two-dimensional cuts 
through parameter space. Note that we used a fiducial survey size of 1 deg^. The cosmological parameters that are not given on 
the axes are evaluated at their fiducial values. In each panel the cross indicates the fiducial set of parameters. The results for the 
data vector D, i.e. the set of tomography power spectra before nulling, are given as shaded contours, where the dark-gray area 
contains 60 % and the light-gray area 80 % of the posterior likelihood. The contour lines indicate the corresponding areas after 
shear-ellipticity nulling, using the full set, i.e. the maximum number of uncontaminated components in D' . Thick lines correspond 
to the 60 % and thin lines to the 80 % credible region. 



the Gram-Schmidt algorithm. This way the vectors correspond- 
ing to the weight functions still fulfill the constraint equation 

To calculate credible regions, the likelihood function in pa- 
rameter space has to be evaluated, which reads under the as- 
sumption of a Gaussian probability distribution function 



and likewise for D' , where Df stands for the data vector, as ob- 
tained for the fiducial model, and p again denotes the set of var- 
ied cosmological parameters. The covariance matrices are only 
evaluated at the fiducial cosmology, as well as the set of weight 



functions ^j^'jW entering the measures in D' . We assume flat 



L{D\p) = \ ^-k[Dip)-D,Yc,-^[Dip)-D,] (g3) 

{2n)-T- VdetO; 
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Fig. 8. Contours of posterior likelihood before and after shear-elUpticity nulling. Shown are the contours for all combinations of 
cosmological parameters out of the set (□„,, Qa, erg, /iioo), the remaining two parameters being marginalized over. The coding of 
areas and curves is the same as in Fig.|7l In addition, the dotted curves enclose the credible regions resulting from using only 
the single optimized weight function as determined in Sect. 13.41 i.e. power spectra of the form n|j*j(^), in D' . As before, thick 
lines correspond to the 60 % and thin lines to the 80 % credible region, while in each panel the cross indicates the fiducial pair of 
parameters. The lower right panel is a detail of the upper left diagram, as outlined by the dashed box. As can be seen here, solid and 
dotted curves nearly coincide. 



priors on the whole range of parameters considered, leading to a 
posterior likelihood 



UD\p) 



(64) 



The boundaries of the four-dimensional grid in parameter space 
over which the sum in the equation above runs and the fiducial 
cosmological parameters are summarized in Table[T] Otherwise 



the setup described in Sect. |3.4| is kept, except for a number of 
adjustments owing to the restrictions in computational power. 
The power spectra are now calculated for Nc = 30 bins in the 
range between ^ = 50 and £ - 10*, making use of the fit for- 
mula f or the non-linear structure evolution bv lPeacock & Doddsl 
(Il996h . In addition, the number of redshift bins is reduced to 
= 20, still ranging from z = to z = 4. 
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Fig. 9. Credible regions in the On, - crs plane for varying redshift binning before and after shear-ellipticity nulling, using the full 
set of measures. In the upper panels the data vectors as given in ( l60b are used; for the results shown in the lower panels the cross- 
correlation power spectra of adjacent bins were also incorporated into both D and D' . The data sets used are otherwise identical, 
except for the number of redshift bins, which is A^; =5 in the left panels, A^; = 10 in the central panels, and A^; = 20 in the 
right panels. The coding of areas and curves is the same as in Fig. [7] Again, the crosses indicate the fiducial set of cosmological 
parameters. 



Table 1. Set of cosmological parameters used for the analysis. 





Pf 


Pmm 




P\ 


PC 




0.3 


0.10 


0.70 


0.20 


0.40 


Q.„ 


0.7 


0.40 


1.00 


0.58 


0.82 




0.9 


0.40 


1.40 


0.80 


1.00 


A 100 


0.7 


0.30 


1.05 


0.60 


0.80 



" Note: values of the fiducial model (pf), limiting values pn,i„ and 
Praax of the parameter plane considered, and lower {pi) and upper limits 
(pa) of the prior applied in the marginalization. 



If all entries in D' of the form n|j^^. are removed, 
only components free of shear-ellipticity correlations remain. 
The loss of information caused by this removal is illustrated in 
Fig.|7] where two-dimensional cuts through the credible regions 
in parameter space, resulting before and after nulling, are given. 
The parameters not shown are evaluated at their fiducial values, 
so that the cross in each panel marks the fiducial model, in this 
case coinciding with the point of maximum likelihood. In the 
^^m - CTg plane, one recognizes the typical banana shape, while 
the Hubble parameter and the density of dark energy are only 
poorly constrained by our setup. As expected, the contours after 
the application of the nulling technique have widened through- 
out. The inner contour line remains outside the light-gray area. 



implying that the probability that a range of parameters contains 
the true cosmological model decreases from 80 % to less than 
60 % after the removal of contamination by intrinsic alignment. 
It is interesting to note that the ratio of the;^'^, i.e. the argument 
of the exponential in ( |63] ), before and after nulling is roughly 
constant over the whole range of parameters considered. 

In Fig.[8]the same set of credible regions is plotted, but here 
the hidden parameters have been marginalized over with flat pri- 
ors within the range indicated in Table[T] The appearance of the 
contours in the Qm - cr^ plane remains similar due to the small 
influence of /zioo and Qa. However, the characteristic banana- 
like shape causes the area of maximum likelihood to be shifted 
away from the fiducial model in panels where one of the param- 
eters Qm or erg is marginalized over Consequently, the poste- 
rior likelihood peaks at lower values of erg in the top center and 
right panel, while Lpost obtains its maximum at lower values than 
Q.m - 0.3 in the central left and lower left panel. Besides that, 
the nulling implies that the degeneracy in /zioo and Qa even in- 
creases, see Fig.|7] which leads to a significant elongation of the 
credible regions along these parameters in the panels of Fig. [8] 
mentioned above. Although this stretching causes lines of equal 
likelihood to even intersect, thereby apparently improving pa- 
rameter constraints in one dimension, it is expected that the total 
area of a credible region increases by shear-ellipticity nulling, 
which corresponds to an overall decrease in constraints. This is 
indeed the case as will be shown below. 
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In addition, we have computed the posterior likelihood for 
nulling with first-order measures alone, i.e. for a vector D' with 
all components apart from those with a subscript [1] removed. 
The resulting contours are also presented in Fig. [8] Within the 
resolution of the graphics and the grid we employed to cover 
the parameter space, the contours for the first-order nulling co- 
incide with the results for the full sefl Only for the more con- 
centrated likelihood in the Q.^ - cr% plane are the first-order con- 
tours located distinctly farther outside, as can be seen in the inlet 
of Fig.|8] The values where the contours are drawn for the two 
setups deviate by less than 3 % in the latter case; for all other 
parameter planes, the deviation of the contour values is approxi- 
mately 1 %. 

We quantify the widening of contours in terms of the 
quadrupole moments of the likelihood function, employing q- 
values as introduced by Kilbinger & Schneider (20()4j). For any 
two-dimensional Ukelihood analysis, one defines 

= ^ L^oAD\p) [pti - Pf,//) [pv - Pf,v) (65) 
P 

for ii,v - 1,2. Then the quantity 

q = Vditfi = ^fiii Q22-Q\2 (66) 

scales with the area of the credible region. Hence, an increase in 
^-value corresponds to a degradation of parameter constraints. 
We compute the relative change in q due to nulling, using both 
the full set and only first-order measures. The results, listed in 
Table|2] confirm the only marginally weaker performance of the 
first-order-only configuration. Depending on the combination of 
parameters, the increase in q ranges between 20 % and about 
50%. 

Table 2. Increase in q'-values for the marginalized credible re- 
gions shown in Fig. [8] 



Pi 


P2 


Ag,o, [%] 






o-% 


29.2 


30.3 




hm 


45.6 


48.3 




£1a 


50.3 


52.5 




o-g 


36.4 


37.5 


fiA 


hm 


21.2 


21.7 




o-g 


26.4 


27.3 



" The relative change in q is given for nulling with the full set (A^,ot) 
and with first-order measures only (Agi). 



Table 3. Increase in ^-values for credible regions, resulting from 
different redshift binning of the survey. 





A^ [%] 


^q' [%z 


5 


192.5 


35.5 


10 


92.7 


17.2 


20 


28.9 


15.6 



" The values Aq correspond to the regions shown in the upper panels, 
the values Aq' to the regions in the lower panels of Fig.|9] 



' Note that 'full set' means the use of the maximum number of new 
power spectra that are not contaminated. 



To assess the importance of a large number of redshift bins, 
we considered a further setup with only two varied parameters 
(Q.ni,o'&), but otherwise identical with respect to the foregoing 
implementation. We used 5, 10, and 20 redshift bins to cover 
the range between z - and z - 4-. The resulting credible re- 
gions before and after nulling, using the full set of measures, are 
given in Fig.|9l upper panels. Concerning the original set of to- 
mography power spectra, the increase in the number of redshift 
bins does not improve the param eter constr aints appreciably, as 
already discussed for instance bv lHij (11999 ). In contrast to this, 
the loss of information due to the removal of potentially contam- 
inated components of the data vector is dramatic for N, = 5 and 
still considerable for A^; = 10; see Table[3]for the corresponding 
changes in ^-values. In the lower panels of Fig.|9]the credible re- 
gions, resulting from the inclusion of the cross-correlation power 
spectra of adjacent bins into the data vectors ( |60] l, are shown (see 
also TableO. The contours tighten substantially, in particular in 
the case A^. = 5, where the probability of a contamination by in- 
trinsic alignment of cross-correlations between neighboring bins 
is low anyway due to the large bin size. With this setup, a num- 
ber of 10 redshift bins already ensures that parameters are still 
well-constrained after nulling. 

5. Conclusions 

In this paper we have elaborated on a purely geometrical method 
to eliminate shear-ellipticity correlations, which constitute a 
source of severe contamination of the cosmic shear signal. Using 
a nulling technique, new observables free of this contamination 
are constructed by suitably weighting shear tomography power 
spectra. The weighting is determined such that the contribution 
to the cosmic shear signal from matter structures potentially 
causing shear-ellipticity correlations is removed by means of the 
characteristic dependence of these coiTelations on redshift. 

Three approaches to obtaining weight functions were inves- 
tigated, which in addition optimize the information content of 
the new observables in terms of the Fisher matrix. The results 
for both analytical and numerical methods are in good agree- 
ment, also for higher orders that are constructed by imposing 
a suitably defined condition of orthogonality. Most notably, the 
analytical ansatz, being computationally simple, is compatible 
in spite of several simplifications. 

Using a set of tomographic power spectra with 20 redshift 
bins, we computed credible regions in a four-dimensional pa- 
rameter space to assess the loss of information due to shear- 
ellipticity nulling. The contours widen significantly with an in- 
crease in ^-values of up to 50 % although reasonably strin- 
gent constraints on cosmological parameters are still possible. 
Besides, it was demonstrated that the use of a smaller sub- 
set of power spectra, excluding those constructed by means of 
higher-order weight functions, yields a practically identical per- 
formance, compared to the full set. 

However, data with a smaller number of redshift bins could 
be shown to almost completely lose its ability to constrain cos- 
mological parameters when intrinsic alignment is taken into ac- 
count (see also BK07), unless the contamination by intrinsic el- 
lipticity correlations of power spectra cross-correlating adjacent 
bins is safely under control. Hence, a large number of redshift 
bins is desirable for cosmic shear studies, in spite of the fact 
that, without the eff'ects of intrinsic galaxy alignment, constraints 
on cosmological parameters are not appreciably improved once 

> 5. These results underline once more the need for both pre- 
cise and detailed redshift information to control systematics in 
cosmic shear. 
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We emphasize that, using the full set of new power spectra, a 
maximization of the information content in the weight functions 
is not necessary. In particular, this means that shear-ellipticity 
nulling as such does not rely on the Fisher matrix formalism 
and the determination of optimal weight functions, nor do the 
credible regions that result from applying the nulling. If only the 
first-order power spectra are employed, parameter constraints do 
depend on the optimization of weight functions. However, the 
first-order results agree well for the three considered approaches, 
so that the change of the likelihood contours due to the use of 
these different methods is expected to be only marginal. 

The moderate change in ^-values indicates that this method 
is in principle suited to inferring cosmological parameters with 
fair precision. It is currently the only truly secure means of 
safely eliminating shear-ellipticity - and in addition, by con- 
struction, intrinsic ellipticity - correlations since no assumptions 
are made about the still uncertain models of these systematic ef- 
fects. An improvement in the performance of shear-ellipticity 
nulling could be achieved by taking advantage of the intrinsic 
ellipticity correlations not being only restricted to galaxy pairs 
which are close in redshift, but also close on the sky, thereby 
allowing for keeping part of the signal from correlations of ad- 
jacent redshift bins. 

As soon as reliable data exists for modeling the intrin- 
sic alignment of galaxies, an intermediate approach between 
strict shear-ellipticity nulling and the full reliance on the shear- 
ellipticity power spectrum could be developed. Since the tidal 
forces acting on a galaxy are caused by the surrounding matter- 
density distribution, coiTelations between matter density and in- 
trinsic ellipticity can be used to determine the tidal field and, 
consequently, the expectation value of the orientation of the in- 
trinsic ellipticity. This additional information, appropriately in- 
corporated into the nulling technique, may improve parameter 
constraints, w hile keeping the influence of model uncertainties 
at a low level. iMandelbaum et alj (1200 6*) have analyzed density- 
ellipticity correlations in the Sloan Digital Sky Survey, using 
galaxies as tracers for the mass; see references therein for fur- 
ther observations. 

A central aspect to be considered before applying the nulling 
technique to real data is the influence of photometric redshift un- 
certain ty since this will fu rther deteriorate constraints (see e.g. 
BK07, lAbdalla et alj|2007l) and possibly enforce a modification 
of the method. Catastrophic outliers in the determination of pho- 
tometric redshifts may cause a leakage of shear-ellipticity cor- 
relations into the measures constructed by the nulling technique. 
Currently, we are not able to judge how severely this leakage will 
affect our method. In a follow-up study, we will investigate the 
demands on photometric redshift accuracy due to nulling in de- 
tail. However, the probability distribution of galaxy distances en- 
ters our method only via the shape noise term in the power spec- 
trum covariance matrix, see STJl . which is merely used for deter- 
mining optimized weight functions and not needed for construct- 
ing the full set of measures, suggesting that shear-ellipticity 
nulling could be fairly robust against incomplete knowledge of 
the distance distribution. Note that, even if very accurate pho- 
tometry were available, the determination of distance in terms 
of redshift would still be limited by the peculiar velocities of 
the observed galaxies and therefore provide an upper limit to the 
number of usable redshift bins. 

Furthermore, it would be of interest to determine the 
decrease in accuracy due to nulling, when determining the 
equation-of-state parameters vvq and w„ of dark energy, because 
upcoming observations will set a focus on these parameters and, 
besides, the results would allow for a direct comparison with 



other methods whose performance is mainly evaluated in terms 
of the dark energy figure of merit. In addition, one may think of 
generalizations such as the application to higher-order statisti- 
cal measures of cosmic shear As far as third-order statistics are 
concerned, neither shear-ellipticity nor intrinsic ellipticity coiTe- 
lations have been subjec t to much investigation so far. Recently, 
ISemboloni et al.l ([2008') have used simulations to demonstrate 
that the three-point shear-ellipticity correlation terms can reach 
up to 10 % of the amplitude of the pure shear signal in a survey 
with a median redshift of 0.7. 

Irrespective of the future practical relevance of the approach 
presented in this study or its possible modifications, the limited 
abilities of this geometrical approach to put constraints on cos- 
mological parameters prove the necessity of improving the un- 
derstanding of galaxy formation and evolution in order to de- 
velop accurate models of intrinsic galaxy alignment. Note that, 
if a good model were known someday, it would be possible to 
even tighten constraints on cosmological parameters since shear- 
ellipticity correlations are another cosmological probe comple- 
mentary to the pure lensing signal. The detailed knowledge of in- 
trinsic alignment would not only allow for the efficient removal 
of contaminations of the cosmic shear signal, but in addition pro- 
vide an interesting means to study the interaction between galax- 
ies and their environment. 
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